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ABSTRACT 


This report presents a detailed description of a computer 
program based on the Range Pair Cycle Counting Method, as given 
in Reference 3, The Range Pair Cycle Counting Method is a procedure 
for generating an analysis spectrum from a given load spectrum. 
Examples are presented where the resulting analysis spectrum will 
be used as input to a crack growth analysis program. 
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SECTION I 


INTRODUCTION 


In crack propagation analysis it is necessary to have a 
correct representation of the load spectrum. A load spectrum 
obtained from tests may not be directly applicable to analysis. 


_The Range Pair Cycle Counting Method is a means of determining 


an accurate analysis spectrum from the actual load spectrum. 

This method is briefly discussed and compared with other counting 
methods in References 1 and 2. A more comprehensive discussion 
which forms the basis for the development of the present computer 
program is given in Reference 3. | 


The computer program treats a load spectrum § as a collection 
of n peaks and valleys designated by xj, i= 1,...,2n, such that if 
x, is a peak then x54, is a valley, 14j<2n-l1. The analysis spectrum 
iz represented by a collection of m cycles {(a,b),},i = 1,...,m, 
such that a; and by are elements of S. The Range Pair Cycle Counting 
Method considers four points (x,,x2,X3,%4) at a time and the 
conditions for counting a cycle (x9,x3) are as follows: 


if Xo > Xy> then a cycle is counted if 
x9 2X4 and X3 2Xy- 
Conversely, if x9< x 1, then a cycle is counted if 
Xo SXy and X35%- 
This method is illustrated in Figure l. 


Thus, starting at the beginning of the load spectrum the first 
four points x}, X29, X3 and x4 are considered. If x9 and x3 meet the 
above conditions, a cycle is defined and these two points are deleted 
from the spectrum. Consequently x4, becomes x2 and the next two 
points of the spectrum are added to again give four points. Counting 
continues until the four points considered do not define a cycle. 
Then x, is omitted from consideration and becomes an element of a 
residue spectrum, The three remaining points are updated, i.e. x9 
becomes Xj, X, becomes x2, x4 becomes X3, and x4 is added sequentially 
from the load spectrum. This process continues until there are only 
two or three points remaining. These points are added to the residue 
spectrum, which is then analyzed in the same manner as the original 
load spectrum. Continuing in this manner a residue spectrum is | 


finally generated which will yield no cycles by the Range Pair 

Cycle Counting Method, This residue spectrum diverges to a maximum 
range and then converges as shown in Figure 2. Cycles are generated 
from the final residue spectrum as follows: Pair the highest peak 
with the lowest valley to form a cycle. Then moving away from this 
cycle in both directions, each successive peak and valley are paired 
together. If there is an extra peak or valley left on either side, it 


-4s omitted, This counting method is illustrated in Figure 2. 


In summary, an original load spectrum is analyzed using the. 
Range Pair Cycle Counting Method to produce an analysis spectrum plus 
a final residue spectrum. This final residue spectrum is then 
analyzed by a pairing technique to determine the remaining cycles, 
which are then added to those previously counted. The result is a 
complete analysis spectrum for use in analytical predictions. 


Fa 


QQ 


SECTION If 


PROGRAM ORGANIZATION 


The Range Pair Cycle Counting program, RPCM, assumes that the 
input load spectrum, S, is defined by n peaks and valleys, (x, i234) and n 
counters k;, i=1,...,n, where k; is a count of the number of eines 
the ith peak and valley are to be repeated sequentially, The program 
then assigns a step number 5] j=1,...2,n to each peak and valley of Se 
Since the analysis spectrum is generated in disjoint parts, i.e. 
from the input load spectrum, from each residue spectrum, and from 
the final residue spectrum, the step numbers are used to sort the 
analysis spectrum relative to the sequencing of the initial load 
Spectrum. Sequence becomes important particularly in crack growth 
analysis. When the counter k is less than 1, as can occur in a flight 
by flight load spectrum, the peak and valley associated with k is not 
analyzed by the program, but is transferred directly into the analysis 
Spectrum and subsequently sequenced relative to its step number. 


The program RPCM is divided into three parts. Each part is 
described below in a step-by-step manner. 


Part I 


1. The initial load spectrum S is adjusted by removing those 
peaks and valleys whose counter k is less than one. 


2. The initial load spectrum S is further adjusted if for 
some i, the ith peak and valley are equal to the (i+1)th peak and 
valley, by maximizing the counter k;. 


3. The Range Pair Cycle Counting Method is now applied to the 
adjusted load spectrum, S. Program RPCM calls Subroutine DECIDE with 
four elements from S. Subroutine DECIDE determines whether a cycle 
is to be generated or whether x, goes to the residue spectrum. Cycles 
are generated in Subroutine CYCGEN. 


Part II 


1. The Range Pair Cycle Counting Method is applied to the © 
residue spectrum, Program RPCM calls Subroutine DECRES with four 
elements from the residue spectrum. Subroutine DECRES determines — 


whether a cycle is to be generated or whether x7 goes to the next _ 


LS a 


\ residue a Cycles | are generated i in “Subroutine ‘CYCRES, 


et ee a ee Fe ae eee eee ee 
_s 
— —_— er ee 
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- - 


| 2. If the current residue spectrum has less than three points 
or if no additional cycles can be generated by the Range Pair Cycle 
Counting Method, proceed to Part III, otherwise return to Step l. 


Part III 


1. The remaining cycles are generated from the final ‘residue 
Spectrum, ) 


2. The analysis spectrum is sorted relative to the sequencing 
of the input load spectrum. 


Card No. 
(Format) 


1 
(8A10) 


2 
(215) 


3,2+0,NPKS+2 _ 
 (5x,3E10. 3) 


SECTION III 


INPUT INSTRUCTIONS 


Variable Name 


TITLE 


NPKS 


NPUNCH 


SIGMAX(T) 


SIGMIN (TI) 


RNCYC (I) 


Definition 


the load spectrum, § 


Number of peaks or valleys 
in the load spectrun, S 


Punch flag 
NPUNCH = O implies the 
analysis spectrum will be 


An alphanumeric description of 


punched in the input format. 


Ith peak of the load spectrum, § 


Ith valley of the load spectrum, § 


counter ky of the Ith peak and 


valley 


SECTION IV 


TABULAR OUTPUT 


Program RPCM gives the following output: 

1. The Lapit load spectrun, Se 

2. The adjusted load spectrum as discussed in Section It, 
3. ie eticiceta and step numbers of Residue Spectrum 1. 


.4, The elements, the step number and the counter k of the cycles 
generated from the adjusted load spectrum. 


5. The elements and step numbers of Residue Spectrum 2, 


6, Step 4 output is repeated plus any additional cycle information 
generated from Residue Spectrum l. 


hs Steps 5 and 6 are repeated for each residue spectrum until the 
final residue spectrum is generated, 


8. All previous cycle output plus any additional cycle information 
generated from the final residue spectrum, 


9, The Range Pair Cycle Counted spectrum, i.e., the analysis 
spectrum, 
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SAMPLE PROBLEMS 


1. 14 Mission C5-A A Spectrum 


A partial plot of the input load spectrum $ is given in : 
Figure 3. The full spectrum is listed on P . It has been ob- 
served that the spectrum listing may not be a good representation 
of the load spectrum since some of the peaks or valley values given 
in the spectrum listing do not match the actual peaks and valleys on 
the load spectrum. This can be illustrated by steps 6 through 9 of 
the spectrum listing. | 


6 8215.0 0.0 1 
7 9146.0 6846.0 b) 
8 6065.0 0.0 — 12 
9 11790.0 8790.0 50 


The load spectrum that these 4 steps would produce is given in 
Figure 4. Now considering the actual peaks and valleys shown in 
Figure 4, steps 6 through 9 should become 


8215.0 0.0 1 
9146.0 6846 .0 4 
9146.0 0.0 1 
6065.0 0.0 11 

0.0 50 


11790.0 879 


These five steps may now be range pair counted according to the 
rules given. The x's on Figure 3 indicate additional places where 
the above type of behaviour occurs. 


The program as written can handle such discrepancies in the 
initial load spectrum. 
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